ARPES spectra of the stripe phase in the 2D t-J model 
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The 2D t-J model with and without t' and t" hopping-terms is studied by exact diagonalization 
on a 5 X 4 cluster, which realizes a hole stripe in y-direction in a spin-Peierls. Next nearest hopping 
terms with a sign appropriate for hole-doped cuprates enhance the stripe formation. The dispersion 
of the quasiparticle-peaks in the single-particle spectrum is in good agreement with bond operator 
theory for hole motion in the spin-Peierls phase, particularly so for realistic values of t' and t" . 
The resulting spectral weight distribution and Fermi surface agree well with experimental ARPES 
spectra on Lai.28Ndo.6Sro.i2Cu04. 

PACS numbers: 71.10.Fd,71.10.Hf,75.40.Mg 
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r~¥ tLe prediction of charged stripes Q and their sub- 
se||)ient experimental verification[2 stand out as one 
oi^he rare instances, where a nontrivial theoretical 
pi^iction for cuprate superconductors was found con- 
s i^t^ nt with experiment. Accordingly, there is cur- 
reiRjly considerable interest in the mechanism leading to 
.ttjen formation and the physical implications of stripes 
■ IIS H II II II IS 113 ■ On the other hand few exper- 
.liAeiital techniques provide such direct experimental in- 
'siHt into the electronic structure of a given compound 
"as^ngle resolved photoemission spectroscopy (ARPES). 
Itjsjtherefore quite natural to look for the fingerprints of 
stripes in ARPES spectra and indeed the results of Zhou 
,' on Lai.28Ndo.6Sro.i2Cu04[ll| are widely consid- 
er^ as strong evidence for stripes. It is the purpose of 
present manuscript to present single particle spectra 
ined by 'computer spectroscopy' on the stripe phase 
qj^e 2D t-J model, presumably the simplest theoretical 
(ifig^ ription of the Cu02 planes in cuprate superconduc- 
t6f^. As will be seen below these result combined with a 
ively crude theory for hole motion in a spin-Peierls 
kground' already give a quite satisfactory description 
ol>«ost of the experimental results. 
Tm t-J model reads 
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denotes summation over pairs of nearest 
sites, Ci^a = Ci,o-(l — rii^a-) and Si and rii de- 
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T^rpby 
n^Hbor 

rt^^ .the operators of electron spin and density at site 
WSpectively. We denote the hopping matrix elements 
ti.^between (l,0)-like neighbors by t, between (1, l)-like 
neighbors by t' and between (2,0)-like neighbors t" , all 
other tij are zero. Throughout we will assume that 
t' /t" = —2, as would be appropriate if the physical ori- 
gin of these terms is hopping via the apex oxgen 2pzS 
orbitalJJJ, and t' /t < 0, as is the case for hole-doped 
compounds. 

The method we apply to study this model is exact di- 
agonalization of finite clusters by means of the Lanczos 
algorithmflj. In a preceding paper^^l we have shown 
that by changing the geometry of the cluster from the 



standard tilted square form to a rectangular one (more 
precisley: to a 5 x 4 cluster) a ground state with a pro- 
nounced stripe like charge inhomogeneity emerges. Here 
we want to discuss the single particle spectra of this state. 
A question to worry about first is, whether the stripes 
survive the additional hopping terms oc t',t". Intuitively 
this is not what one would expect, because additional 
hopping terms increase the mobility of the holes, whence 
any spatial inhomogeneity should be washed out more 
efficiently. Surprisingly enough, the numerics show, that 
exactly the opposite is happening: the additional hop- 
ping terms even slightly enhance the charge inhomo- 
geneity. This is demonstrated in Table ^ which com- 
pares the static density correlation function g]j(R) = 
'^j{njnj-\--R.) for vanishing and finite t' and t" . Next 
nearest-neighbor hopping terms with the proper sign for 
hole-doped cuprates thus seem to have a stabilizing effect 
on stripes - if any. 

Next we adress a special feature of the 5x4 cluster, 
which will be essential to understand the hole dynamics, 
namely the presence of spin-Peierls dimerization even at 
half filling. TablelTTlshows that the static spin correlation 
function is strongly anisotropic, with singlet-bonds pre- 
dominantly in y-direction. Since the boundary conditions 
in 5 X 4 frustrate the Neel order, they apparently stabi- 
lize the energetically close spin-Peierls phase. In fact, 
the GS energy of the 5x4 cluster is only marginally 
higher than that of the square-shaped x cluster 
(— 1.165J/site vs. — 1.191J/site). Clearly, this is a con- 
firmation of the proposal by Read and SachdevfTEj that 
the transition to a spin-Peierls phase is a likely instability 
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0.298 0.139 0.035 
0.239 0.144 0.041 
2.000 0.049 0.018 


0.303 0.151 0.037 
0.261 0.144 0.043 
2.000 0.016 0.010 
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TABLE I: Static density correlation function (;_d(R), 5x4 
cluster with 2 holes, J/t = 0.5. Other parameters are t' — 
t" = (left panel) and t' /t = -0.4, t" /t = 0.2 (right panel). 
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FIG. 1: Left; Photoemission spectrum (PES) for the half- 
filled 5x4 cluster, J/t = 0.5. , t' /t = -0.2, t" /t = 0.1. 
Right: Dispersion of the quasiparticle peak as extracted from 
the numerical spectra (top) compared to the theoretical single 
hole dispersion ei(k) in the spin-Peierls phase (bottom). 



of the S = ^ 2D Heisenberg antiferromagnet. 

Clear evidence for the spin-Peierls nature of the half- 
filled ground state can be seen in the single particle spec- 
tral function ^(k, ui) (see Ref. |l3j for a definiton), which 
is shown in Figure ^ It is immediately obvious that this 
differs markedly from the familiar dispersion for a hole 
in an antiferromagnet: whereas for hole motion in a Neel 
state the top of the ARPES spectrum is at (|, |)[ll,|T3|, 
the dispersion seen in the half-filled 5x4 cluster has its 
maximum at (^,7r/2) - which probably means (tt, 7r/2) 
in the infinite system. Another notable feature is the 
symmetry of the dispersion under the exchange {k^, 0) 
{kx,Tr) - which is exactly what one would expect from 
the doubling of the unit cell by spin-Peierls order with 
dimers in y-direction. To be more quantitative, let us dis- 
cuss the single-hole dispersion in the spin-Peierls phase. 
Starting from a product state of singlets, which cover the 
bonds of the lattice in the form of a columnar pattern, we 
derive a Hamiltonian for the motion of singly occupied 
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0.222 -0.162 0.052 
-0.389 0.173 -0.059 
0.750 -0.276 0.061 
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TABLE II: Static spin correlation function for the half- filled 
5x4 cluster. 



dimers. The singlet state on a bond formed by the sites 
(1,2) is |s) = -^(c| I — c[ |)|0) . A dimer with a 
single hole can be in either the bonding or antibonding 
state: = ^'W*^) (gt^^ j- c\^^)\0). Introducing the 

'creation operator' h^. ^. = |±,a)(s| we can - by straight- 
forward generalization of Refs. [la . ITgl l20l| - write down 
the following Hamiltonian describing the motion of these 
effective Fermions: 

H = 5]6+(k)/i1;,,>+^k,. + e-(k)/i^k>-,k,. 

k,(T 

+ (nk) hl ^^h_ ,,^, + H.c)j 

6±(k) Ti + t(^cos(fc.)±£^5|M^ 

±t' cos{k^) (1 + cos(2fcy)) 

+t" (cos(2fc:r) + cos(2fcy)) 

it 

y(k) = --sm{2ky) ~ it' cos{k^)sm{2ky). (1) 

Here the coupling of the dimer- Fermions to triplet- 
excitations of the bonds flfll l2(ij has been neglected 
for simplicity. Diagonalizing the Hamiltonian we 
obtain the dispersion relation eQ,(k) and the quasipar- 

ticles 7^ = UQ,k/i+,q,cr + ^'a.k/i-,q,o-> whcrc a G {1, 2}. 
Figure Q shows that there is good agreement between 
the numerical peak dispersion in the 5x4 cluster and 
our simple theory. The main differences are the flatten- 
ing of the cluster dispersion near the band maximum at 
(^,7r/2) and the smaller bandwidth in the numerical 
spectra, which is probably due to the coupling of the 
Fermions to triplet excitations [l^ ITqI l20l|. Taking into 
account the simplicity of our calculation, however, the 
agreement is quite satisfactory. Summarizing our discus- 
sion so far we have seen that the 5x4 cluster shows 
a single-hole dispersion that differs markedly from that 
for a Neel background, but is in good agreement with 
a simple bond-operator calculation for hole-motion in a 
spin-Peierls phase. This reflects the fact that the frus- 
tration of Neel order along the odd-numbered side of the 
5x4 has driven the system into the spin-Pcicrls phase of 
the 2D Heisenberg antiferromagnet, consistent with ar- 
guments given by Read and Sachdev[l5l|. 
While the dimerization clearly renders the 5x4 clus- 
ter unsuitable to describe undoped compounds such as 
Sr2CuO2Cl20|, it makes the interpretation of the spec- 
tra for the striped ground state at finite doping a lot 
easier - as will be seen now. Ignoring for the moment the 
formation of a hole stripe as well as the fact that the 
dimer- Fermions ^± q o- actually obey a hard-core con- 
straint, one would expect, that the doped holes accu- 
mulate near the top of the single-hole dispersion, thus 
forming to simplest approximation a cigar-shaped hole 
pocketfl^ l20| centered at (^, -1) (or rather (tt, |) in the 
thermodynamical limit). Figure [3 which shows the sin- 
gle particle spectrum for the two-hole GS for different 
t' and t" , demonstrates that this is indeed exactly what 
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FIG. 2; Single particle spectrum for the 5x4 cluster with 
two holes. Parameter values are J/t = 0.5. , t' — t" — (Left) 
and J/t = 0.5. , t' /t = -0.2, t" /t = O.f (Right). Only the 
parts near the Fermi energy Ep/t ~ 1.5 are shown. 



FIG. 3: Left: Single-particle spectral function from bond- 
operator theory for t' /t — —0.1, t" /t = 0.05, the Fermi energy 
(verical line) corresponds to a hole density of 0.16. 
Right: Single hole dispersion (top) and momentum distribu- 
tion n(k) (bottom) for 5x4 and the same parameters. 



happens. There is a clear analogy between the spectra 
for the doped case and half-filling, the dispersion of the 
quasiparticle peak being essentially unchanged. There 
are two major changes, both of whom could have been 
expected: first, peaks with a higher binding energy be- 
come rather diffuse, which is nothing but the familiar 
Landau-damping. Second, the peaks at (-^if-) and (to 
a lesser extent) at f ) cross to the IPES spectrum. 
It should be noted that due to finite-size effects there is 
always a significant gap between the PES and IPES spec- 
trum of finite clusters - after all the electron numbers of 
initial and final state differ by a finite fraction (10% in 
the present case). It is therefore impossible to decide, 
whether the gap is simply a finite-size effect or due to 
the stripe formation. The only deviation from this ideal 
rigid-band behaviour is the appearance of high-energy 
IPES peaks along (0,7r) (^,7r). The interpretation 
of these peaks, however, is straightforward: in inverse 
photoemission, an electron is necessarily inserted into a 
dimer occupied by a single electron. The spins of the 
two electrons then can couple either to a singlet - which 
means the IPES process leads back to the spin-Peierls 
'vacuum' - or to a triplet - which means the IPES process 
leaves the system in a spin-excited state. The IPES peaks 
along (0, 7r) (tt, tt) presumably originate from the lat- 
ter process (we note that exactly the same holds true also 
for the 'usual' ground state of the t-J model, see 0|). 
Taken together, the data presented so far demonstrate, 
that the spin-Peierls order in the striped phase is the 



key to understand its single particle spectra. Neglecting 
subtleties such as the possible formation of a Luttinger 
liquid along the stripes 0, the possible condensation of 
d-like hole pairs along the stripes 1231 or the formation of 
various kinds of order parameter[20j, the spectra show, 
that the system can be described by the dispersion for 
a single hole in a spin-Peierls background being filled up 
with holes. Thereby the momentum (tt, ^) where the 
band crosses the Fermi energy is independent of t' and 
t" - the reason is simply that irrespective of t' and t" 
the dispersion is symmetric under reflection by the line 
(0, ^) (tt, ^) due to the underlying spin-Peierls order. 

We now want to use this (rather oversimplified) sce- 
nario to discuss the experimental ARPES spectra on 
Lai.28Ndo.6Sro.i2Cu04llll. To compute the full ARPES 
spectrum we need the following representation of the elec- 
tron annihilation operator: 

Ck,. = (cos(^) hi^_^^, - .sin(^) . (2) 

Thereby q is the 'backfolded version' of k, so as to take 
into account that the Brillouin zone of the spin-Peierls 
phase is [— 7r,7r] x [— ^, Equation (jSJ is readily veri- 
fied by taking matrix elments of both sides between \s) 
and |±,cr). The spectral weight of a given quasiparticle 
branch then is w = |u* ^ cos(^) + iw* sin(^)p. In 
FigureOlthe spectra for the symmetric momenta {k^, ky) 
and {ky,kx) have been averaged, as would be appropri- 
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ate for an ARPES experiment on a compound with do- 
mains of different singlet direction|llj. The experimental 
ARPES spectrum should be compared to the parts of the 
spectrum below the Fermi energy, indicated as the ver- 
tical line in Figure |31 Along the (1, 1) direction there is 
a band dispersing upwards and disappearing halfway be- 
tween (0, 0) and (tt, tt). This is actually not a Fermi level 
crossing - the dispersion of the band actually bends down- 
wards again after passing through (f,f) - but rather 
a vanishing of the spectral weight due to a destructive 
interplay between the coefficients (Ma,k,i'a,k) and the 
'form factors' {cos{^) , i sm{^)) of the dimer fermions. 
The situation along (1, 1) thus is similar to the 'remnant 
Fermi surface' in the half- filled compounds p^. Along 
(tt, 0) (tt, tt) there are two 'real' Fermi surface cross- 
ings, which are symmetric around (tt, -1). Again due to 
matrix element effects, the one between (tt, ^) and (tt, tt) 
has small spectral weight, which would probably render 
it unobservable in an ARPES experiment. 'Spectroscop- 
ically' the system thus would look very much like having 
a single Fermi surface sheet near (vr, 0) running roughly 
parallel to (1,0) (see FigureEJ and disappearing as (§, J) 
is approached. Compared to experiment, the band cross- 
ing along (tt, 0) (tt, tt) is too far away from (tt, 0) - it 
should be noted, however, that in the actual stripe phase 
the hole density is inhomogeneous, so that the hole den- 
sity 'within' the stripes - which probably determines the 
area of the cigar - is higher than average. Clearly, this 
would shift the Fermi surface crossing towards (tt, 0). No 
crossing should be seen along (1,1), only a band dis- 
persing upwards and disappearing somehwere below Ep . 
The agreement with experiment would be almost perfect 
if the band portion near (tt, 0) were slightly moved up- 
ward compared to the band maximum near (§,§)- such 
fine details are probably beyond our very simple bond- 
operator calculation. To compare to the data of Zhou 
et al. Figure |21 also shows the momentum distribution 



n(k) = A-{]<.,LLj)du! obtained from the cluster diag- 
onalization. This shows some anisotropy, with in par- 
ticular a 'ridge' of almost constant n(k) running along 
(1,0). The cross-shaped character is not as pronounced 
as in the experiment, but this is simply due to the 'par- 
tial occupation' of the momentum (^, §), which crosses 
only partially from PES to IPES, see Figure |21 nd thus 
retains a relatively large n(k). Nevertheless one may say, 
that when combined with the simple bond-operator the- 
ory, the numerical data give a good description of the 
experimental data on Lai.28Ndo.6Sro.i2Cu04. 
In summary, numerically exact diagonalization results 
show, that the 2D t-J model with next-nearest neigh- 
bor hopping terms appropriate to describe hole doped 
cuprates has a spin-Peierls phase which leads to the for- 
mation of pronounced hole-stripes. Given that stripes 
are an experimental fact in cuprates and that the den- 
sity matrix renormalization group calculations on much 
larger clusters of the t-J model 6, 9J have also shown clear 
evidence for stripes, one may be confident that these hole- 
stripes are indeed representative for those in the bulk 
system. The single particle spectrum in the stripe phase 
then is found to be in good agreement with a simple 
bond-operator theory for hole motion in a spin-Peierls 
phase, thus providing further evidence for the intimate 
relationship between spin-Peierls odcring and stripe for- 
mation. Upon doping, holes accumulate near the top 
of the single-hole dispersion, to simplest approximation 
forming cigar-shaped pockets centered on the corner of 
the spin-Peierls Brillouin zone at (tt, -1) - the latter in 
full agreement with bond-operator theorv[l9l l20j. The 
notion of a Fermi surface should not be taken too lit- 
eral, because close to Ep the stripe formation is likely to 
change this simply free-particle picture drastically, but 
all in all the quasiparticle dispersion and spectral weight 
distribution of the stripe phase as seen in the simulations 
are in good agreement with ARPES experiments. 
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